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Abstract 

We study the decay of D° and Df mesons into five-body final states including 
a Kg and report the discovery of the decay mode Df — > KgKg^^ir^. The 

branching ratio for the new mode is ^r^ffc*!"!^ = 0.102±0.029±0.029. We 

also determine the branching ratio of r ^,7* f f g L ) 7r , - !*" - = 0.095±0.005±0.007 as 

well as an upper limit for r poIl^^+ J+ J- J- ) < 0.054 (90% CL). An analysis of 
the resonant substructure for D° — > A"g7r + 7r + 7r _ 7r _ is also performed. 
PACS numbers: 13.25.Pt, 14.40Lb 



More information on multibody final states in the charm sector is an essential 
ingredient for our ability to model decay rates and to further increase our 
understanding of the decay process in heavy quark systems. This is particularly 
important for the decays where a substantial part of its hadronic decay 
rate is still not identified. In this paper we extend our work [1] on four-body 
decays involving a K s to five-body decays involving a Kg. We have already 
published results on all charged five-body modes [2]. The FOCUS collaboration 
presents the first evidence of the decay mode Df — > KgKgTT + 7i + 7i^ , measures 
an inclusive branching ratio for the mode D° — > K^ir + n + n~ n~ relative to 
D° — > Kgir + n~ and places an upper limit on the mode D° — > KgK~-n + n + n^ . 
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Finally we present the first resonant substructure analysis of the decay mode 
D° -> K° s n + n+n-n-. 

The data were collected during the 1996-1997 fixed target run at Fermilab. 
Bremsstrahlung of electrons and photons with an endpoint energy of approx- 
imately 300 GeV produces photons which interact in a segmented beryllium- 
oxide target to produce charmed particles. The average photon energy for 
events which satisfy our trigger is ~ 180 GeV. Charged decay products are 
momentum analyzed by two oppositely polarized dipole magnets. Tracking is 
performed by a system of silicon vertex detectors [3] in the target region and 
by multi-wire proportional chambers downstream of the interaction. Particle 
identification is performed by three threshold Cerenkov counters, two electro- 
magnetic calorimeters, a hadronic calorimeter, and two muon systems. 

Five-body D° and decays are reconstructed using a candidate driven vertex 
algorithm [4]. A decay vertex is formed from the reconstructed charged tracks. 
The Kg is also reconstructed using techniques described elsewhere [5]. The 
momentum information from the Kg and the charged tracks is used to form 
a candidate D momentum vector, which is intersected with other tracks to 
find the production vertex. Events are selected based on several criteria. The 
confidence level for the production vertex and for the charm decay vertex 
must be greater than 1%. The reconstructed mass of the K s must be within 
four standard deviations of the nominal Kg mass. The likelihood for each 
charged particle to be a proton, kaon, pion, or electron based on Cerenkov 
particle identification is used to make additional requirements [6]. For pion 
candidates we require a loose cut that no alternative hypothesis is favored 
over the pion hypothesis by more than 6 units of log-likelihood. In addition, for 
each kaon candidate we require the negative log-likelihood kaon hypothesis, 
Wk = —2 ln(kaon likelihood), to be favored over the corresponding pion 
hypothesis W n by W n — Wk > 2. We also require the distance between the 
primary and secondary vertices divided by its error to be at least 10. Finally, 
in order to reduce background due to secondary interactions of particles from 
the production vertex, we require the secondary vertex to be located outside 
the target material. 

For individual modes we apply additional analysis cuts. Due to the large com- 
binatoric background for D° — > Kg7i + ii + ii^7T~ , we increase the separation 
requirement of the secondary vertex from being just outside the target mate- 
rial to two standard deviations from the edge of the target material. Figure 
1(a) shows the Kg7i + n + n^'K~ invariant mass plot for events that satisfy these 
cuts. The distribution is fitted with a Gaussian for the D° signal (1283±57 
events) with the width and mass floated and a first degree polynomial for 
the background. Figure 1(b) shows the Kgn + n + 7Y^7Y^ invariant mass plot for 
events originating from a D* + — > D°n + decay. 
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Fig. 1. Invariant mass distributions for (a) Kg7r + 7r + 7r~7r~, (b) K^ir + ir + ir~ Tt~ for D* 
tagged events, (c) KgKgTr + ir~TT ± , and (d) KgK~ir + ir~ir + . The fits are described 
in the text. 

The -D+ — > K s K s 7!- + n + n^ mode is difficult to detect due to the relative inef- 
ficiency of Kg reconstruction and that most of the time only the three pions 
define the secondary vertex. The confidence level that a pion track from the 
decay vertex intersects the production vertex must be less than 2%. We also 
require a reconstructed momentum of greater than 25 GeV/c. Figure 1(c) 
shows the KgKg7i + 7i + n^ mass plot for events which satisfy these cuts. This 
is the first observation of this mode. We fit with a Gaussian (37±10 events) 
with mass and width allowed to float and a second degree polynomial for the 
background. 

The decay D° — > K s ) K~'k + 'k + 'k~ is Cabibbo suppressed, and we do not observe 
a signal in this mode. Thus we choose our analysis cuts by maximizing the 
quantity S/ \Z~B~, where S is the fitted yield from our Monte Carlo simulation of 
the mode, and B is the number of background events in the signal region from 
data. Based on this optimization we require a reconstructed D° momentum 
of greater than 50 GeV/c. We also require the D° come from a D* + decay, 
that is 0.142 GeV/c 2 < M D *+ - M D o <0.I49 GeV/c 2 . Figure 1(d) shows the 
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resulting K S K it + -K + it invariant mass plot. As there is no apparent signal 
we report an upper limit branching ratio. 

Table 1 

Branching ratios for modes involving a Kg. All branching ratios are inclusive of 
subresonant modes. 



Decay Mode 


Branching Ratio 


r(D°->K°iv+Tr+-n- w~ ) 

r(i)o_^o 7r + 7r -) 


0.095±0.005±0.007 
0.102±0.029±0.029 
< 0.054 (90% C.L.) 


r(D+^K0K-n+TT+) 
TiD^K^K-n+n+Tr-) 


r(D°^K^7T+7T+n-TT-) 



We measure the branching fraction of the D° — > K s 7i + 7i + 7i n mode rela- 
tive to D° -> K s it + it- ". The relative efficiency is determined by Monte Carlo 
simulation. The K s n + n^ and K s 7i + n + 7Y^n^ channels are produced as an in- 
coherent mixture of subresonant decays based on PDG information [7] and our 
analysis described below, respectively. We measure the — > K s K s 7t + 7t + tt~ 
mode relative to — > K s K~tt + tt + . We test for dependency on cut selection 
in both modes by individually varying each cut. The results are shown in Table 
1, and we compare our measurement of the D° — > K s 7!- + n + n^n^ branching 
ratio with previous measurements in Table 2. 

We studied systematic effects due to uncertainties in the reconstruction effi- 
ciency, in the unknown resonant substructure, and on the fitting procedure. To 
determine the systematic error due to the reconstruction efficiency we follow 
a procedure based on the S-factor method used by the Particle Data Group 
[7]. For each mode we split the data sample into four independent subsamples 
based on D momentum and on the period of time in which the data was col- 
lected. These splits provide a check on the Monte Carlo simulation of charm 
production, of the vertex detector (it changed during the course of the run), 
and on the simulation of the detector stability. We then define the split sam- 
ple variance as the difference between the scaled variance and the statistical 
variance if the former exceeds the latter. The method is described in detail 
in reference [11]. In addition, we split the data sample into three independent 
subsamples based on the location and geometry of the K s decay. We then cal- 
culate the K s reconstruction variance using the same procedure described for 
the split sample variance. We also vary the subresonant states in the Monte 
Carlo and use the variance in the branching ratios as a contribution to the 
systematic error. We also determine the systematic effects based on different 
fitting procedures. The branching ratios are evaluated under various fit condi- 
tions, and the variance of the results is used as an additional systematic error. 
Finally, we evaluate systematic effects from uncertainty in the absolute track- 
ing efficiency of multi-body decays using studies of D° — > K~7t + 7t + tt~ and 
D° — > K~n + decays. The systematic effects are then all added in quadrature 
to obtain the final systematic error. 
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Table 2 

Comparison of this measurement of D° — > i^27r + 7r + 7r _ 7r _ mode to previous mea- 
surements. 



Experiment Events 



r(D"^K°1T + 1T+1T-1T-) 
r(D°~*K i l.TT+TT-) 



E831 (This Measurement) 1283 0.095±0.005±0.007 

PDG Average[7] 0.107±0.029 

ARGUS[8] 11 0.07±0.02±0.01 

CLEO[9] 56 0.149±0.026 

E691[10] 6 0.18±0.07±0.04 

We do not observe a signal in the decay D° — > KgK^n + 7i + 7Y^ and we calculate 
an upper limit for the branching ratio with respect to D° — > Kgn + n + 7Y^7Y^ . 
We evaluate the upper limit using the method of Rolke and Lopez [12]. We 
define the signal region as being within ±2er of the nominal D° mass, and the 
two sideband regions as 4-8(7 above and below the D° mass. We observe 3 
events in the signal region and 6 events in the sidebands, corresponding to an 
upper limit of 5.02 events (@90% CL). 

We study systematic effects for this channel from cut variation and resonant 
substructure, and include these in our determination of the upper limit using 
the method of Cousins and Highland [13]. We determine the systematic error 
from cut variation by individually varying each cut, fitting the resulting dis- 
tribution, and taking the variance between each branching ratio measurement 
as our systematic error. We also study systematic effects from our uncertainty 
in the resonant substructure of the mode by varying the subresonant states 
included in the Monte Carlo simulation, and used the variance in the resulting 
branching ratios as our systematic error. These two systematic effects are then 
added in quadrature to give a final relative systematic error of 26%. 

We then determine the increase in our upper limit based on the equation: 

Atf = ±£/V U + b ~ S 



2 sys U + b 

where U is the original upper limit of events, a sys is the percent systematic 
error determined above, b is the number of events observed in the sideband 
region, and s is the number of signal events observed. We calculate an upper 
limit of 5.64 events, corresponding to an upper limit for the branching ratio 
of: 



T(D° -> K° s ir+n+n-n- 



< 0.054 (@90% CL). 



We have studied the resonance substructure in the decay D° — > Kg7Y + n + n n . 
We use an incoherent binned fit method [14] developed by the E687 Collabo- 
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ration which assumes the final state is an incoherent superposition of subreso- 
nant decay modes containing vector resonances. A coherent analysis would be 
difficult given our limited statistics. For subresonant decay modes we consider 
the lowest mass (K s it~) and (7r + 7r~) resonances, as well as a nonresonant chan- 
nel: K*-ti + ti + ti-, K° s p°7t + iT~, K*-p n + and (K° s it + it + it-it-) nr . All states not 
explicitly considered are assumed to be included in the nonresonant channel. 

For the resonant substructure analysis of D° — > Kgn + n + n^n^ we place ad- 
ditional cuts to enhance the signal to background ratio. We require the confi- 
dence level that a track from the decay vertex intersects the production vertex 
be less than 8%. We also require the D° to come from a D* + decay, that is 
0.144 GeV/c 2 < M D , + - M D o <0.148 GeV/c 2 , in order to reduce background 
and distinguish between D° and D°. Fig. 1(b) shows the Kg7i + n + 7Y^n^ invari- 
ant mass plot for events which satisfy these cuts. We then determine the accep- 
tance corrected yield into each subresonant mode using a weighting technique 
whereby each event is weighted by its kinematic values in three submasses: 
(Kg7r~), (7r + 7r~), and (7r + 7r + ). No resonance in the (tt + 7t + ) submass exists, 
but we include it in order to compute a meaningful x 2 estimate of the fit. 
Eight population bins are constructed depending on whether each of the three 
submasses falls within the expected resonance (In the case of ir + ir + , the bin is 
split into high and low mass regions). For each Monte Carlo simulation the bin 
population, rii, in the eight bins is determined and a matrix, Tj a , is calculated 
between the generated states, a, Monte Carlo yields, Y a , and the eight bins i: 

a 



The elements of the matrix, T, can be summed to give the efficiency for each 
mode, e a \ 

i 

The Monte Carlo determined matrix is inverted to create a new weighting 
matrix which multiplies the bin populations to produce efficiency corrected 
yields. The weight includes the contributions from the four combinations we 
have for each event. Each data event can then be weighted according to its 
values in the submass bins. Once the weighted distributions for each of the 
four modes are generated, we determine the acceptance corrected yield by 
fitting the distributions with a Gaussian signal and a linear background. Using 
incoherent Monte Carlo mixtures of the four subresonant modes we verify that 
our procedure is able to correctly recover the generated mixtures of the four 
modes. 

The results for Kgn + 7c + 7T~7T~ are summarized in Table 3. The four weighted 
histograms with fits are shown in Fig. 2, where Fig. 2(e) is the weighted distri- 
bution for the sum of all subresonant modes. The goodness of fit is evaluated 
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Table 3 

Fractions relative to the inclusive mode for the resonance substructure of the D° — > 
Kgir + ir + Tr~ir~ decay mode. These values are not corrected for unseen decay modes. 

Subresonant Mode Fraction of i"Lg7r + 7r~7r + 7r~ 



{K^TT + TT + TT-TT-) NR < 0.46 @90% CL 

K*-ir + ir + ir- 0.17±0.28±0.02 

K%p°tt + it- 0.40±0.24±0.07 

K*-p°TT+ 0.60±0.21±0.09 
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Fig. 2. i^^7r + 7r + 7r 7r weighted invariant mass for (a) {K^it + it + it it )nr, (b) 
K*~ir + Tr + ir~ , (c) i<L£ / o 7r + 7r~, (d) K*~p°ir + , (e) Inclusive sum of all four modes. 

by calculating a x 2 for the hypothesis of consistency between the model pre- 
dictions and observed data yields in each of the 8 submass bins. The calculated 
X 2 is 9.7 (4 degrees of freedom), with most of the x 2 contribution resulting 
from a poor Monte Carlo simulation of the n + 7i + spectrum in the nonresonant 
channel. 

We observe results similar to previous studies of five-body charm decays, with 
a small nonresonant component and the dominant mode of the form vector- 
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vector-pseudoscalar. Such a result has been predicted by theoretical discussion 
of a vector-dominance model for heavy flavor decays [15], which suggests that 
charm decays are dominated by quasi-two-body decays in which the W ± im- 
mediately hadronizes into a charged pseudoscalar, vector or axial vector me- 
son. Results consistent with the vector-dominance model have already been 
seen by FOCUS in five-body decays [2]. Such theoretical discussion raises the 
possibility that the resonant substructure for the decay D° — > Kgir + 7T + 7T~7T~ 
is dominated by the quasi-two-body decay K*~af. To test this hypothesis 
we generate Monte Carlo simulations of this decay, assuming the af has a 
width of 400 MeV/c 2 and decays entirely as an S-wave to p°ir + , and use our 
subresonant analysis procedure explained above. We observe yield fractions in 
each of the subresonant modes similar to the reported fractions from the data, 
suggesting our results are consistent with the decay being dominated by the 
K*~af subresonant state. 

In conclusion we have measured the relative branching ratios of many-body 
hadronic modes of D° and involving a Kg decay and have presented the 
first evidence of the decay mode Df — > KgKg7!- + n + n^ . We have also per- 
formed an analysis of the resonant substructure of the decay D° — > KsTt + it + it~it 
Finally we have placed an upper limit on the relative branching fraction of 
the Cabibbo suppressed decay D° — > KgK~7r + 7r + 7r~ . 
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Tecnologico, CONACyT-Mexico, and the Korea Research Foundation of the 
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